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Abstract of DE1 98291 70 

The invention relates to a method for computing and displaying 2D- or 3D- 
spatial differences for medical diagnosis, research or medical equipment. 
According to said method, it is possible to compute, store and display 
relative positions of volume elements of at least one object, said elements 
being detected through digitisation achieved by an image formation 
method, whereby at least one reference structure and generally into 
complete volume elements of said reference are formed. The invention is 
characterised in that the volume elements of at least two objects are 
detected by a computer-aided system, said elements resulting from 
digitisation of the image formation method. The paths followed by the 
volume elements of the initial volume relatively to the reference volume 
elements are computed by a linearised displacement model, by a 
multiframe digital method adapted to said displacement model, by forces 
based on grey values and/or with several resolution stages of the digitised 
volumes. 
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Berechnung und Darstellung 2- oder 3-dimensionaler Deformationen raumlicher Strukturen 

Verfahren zur Quantifizierung und Darstellung von 2- wiridiche Bewegnns 

oder3-dimensiona!en Raumunterschieden fur die medizi- 
nische Diagnostik, Forschung oder Medizintechnik, mit 
dem die relativen Positionen der aus der Digitalisierung 
bildgebender Verfahren erfa&ten Vblumenelemente min- 
destens eines Objekts zu mindestens einer Referenzstruk- 
tur und i. A. zu samtlichen Volumenelementen dieser Re- 
ferenz berechnet, gespeichert und da rgestellt werden, da- 
durch gekennzeichnet, daf5 die aus der Digitalisierung 
bildgebender Verfahren resultierenden Vblumenelemen- 
te mindestens zweier Objekte von einem computer-ge- 
stutzten System erfa St werden, wobei die zuruckgelegten 
Wege der Vblumenelemente des Ausgangsvolumens auf 
die entsprechenden Referenz-Volumenelemente durch 
ein linearisiertes Bewegungsmodell, ein auf dieses Bewe- 
gungsmodell angepafctes numerisches Mehrgitter-Ver- 
fahren, grauwertbasierenden Kraften und/oder mit meh- 
reren Auflosungsstufen der digitalisierten Volumina be- 
rechnet werden. 
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Beschreibung 



Das Verfahren dient zur Berechnung und Darstellung von 
2- oder 3-dimensionalen Raumunterschieden (Deformati- 
onsfelder) zur raumlichen Nonnierung von biologischen 5 
Objckten zur Anwendung flir die medizinische Diagnoslik, 
Forschung oder Medizintechnik. Die Gesamtheit der Rau- 
munterschiedc zweier Objekte in eincm gemeinsamen Ko- 
ordinatensystem wird im Folgenden als Deformationsfeld 
bezeichnet. Hierbei werden nicht die Objekte sondern die io 
relauven Positionen der aus der Digitalisierung bildgeben- 
der Verfahren resuluerenden Volumenelemente mindestens 
eines Objekis zu mindestens ciner Refercnzstruktur und 
i. A. zu samtlichen Volumenelementen dieser Referenz be- 
rechnet und dargestellt. Zur Berechnung der 2- oder 3-di- IS 
mensionalen Deformationsfelder verwendet das Verfahren 
ein computergestiitztes System zur Aufnahme der analog zu 
digital gewandelten Signale (Volumenelemente) bildgcben- 
der Verfahren. Nach der Erfassung werden entsprechende 
Volumenelemente im Ausgangs- und Referenzvolumen mit 20 
Hilfe eines Bewegungsmodells analysiert und zugeordnet. 
Die Zuordnung erfolgt dabei durch grauwertbasierende 
Krafte zwischen dem Ausgangs- und Referenzvolumen, wo- 
mit der zuriickgelegte Weg jedes Volumenelementes des 
Ausgangsvolumens mit Hilfc von numerischen Mehrgitter- 25 
Verfahren berechnet werden kann. Der von jedem Volumen- 
element zuriickgelegte Weg wird in einem MeBwertspeicher 
erfaBt. Daran anschlieBend wird die punktweise Verschie- 
bung des Ausgangsvolumens mit Hilfe des Deformations- 
feldes durchgefuhrt und das Ergebnis von einem Speicher 30 
erfaBt oder einer Visualisierungseinheit zugefuhrt. Das De- 
formationsfeld kann auch durch Anwendung an einem 2- 
oder 3-dimensionalen regelmaBigen Gitter auf einer Visuali- 
sierungseinheit sichtbar gemacht werden und dient zur stati- 
stischen Analyse der Variabilitat von raumlichen Strukturen 35 
(interindividuelle Anpassung) oder der Uberlagerung meh- 
rerer Objekte, die mindestens durch zwei verschiedene bild- 
gebende Verfahren dargestellt sind (intra- und interindividu- 
elle Anpassung). Insbesondere das Problem der automati- 
schen, geomelrischen Zuordnung von Volurnina zweier 40 
bildgebender Verfahren zur Informationsverdichtung funk- 
tioneller und morphologisc her Parameter wird mit dem Ver- 
fahren gelost. 
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Die Erfindung dient zur Berechnung und Darstellung von 
2- oder 3-dimensionalen Raumunterschieden (Deformati- 
onsfelder) fur die medizinische Diagnosuk, Forschung oder 
Medizintechnik. Hierbei werden nicht die Objekte sondern 50 
die relativen Positionen der aus der Digitalisierung bildge- 
bender Verfahren rcsultierenden Volumenelemente minde- 
stens eines Objekts zu mindestens einer Referenzstruktur 
und i. A. zu samtlichen Volumenelementen dieser Referenz 
berechnet und dargestellt. Es ist bereits bekannt, bildge- 55 
bende Verfahren in der radiologischen Diagnostik einzuset- 
zen, um normalc und pathologisch vcrandcrte Strukturen 
und Funktionen des Korpers bildmaBig zu erfassen und zu 
analysieren. Die verschiedenartigen Verfahren basieren auf 
einem gemeinsamen Grundprinzip: verschiedene "Informa- 60 
tionstrager" wie z. B. sichtbares Licht, Rontgenstrahlen, Ul- 
traschall, RadionukUde, Magnetfelder treten in Wechselwir- 
kung mit dem Korper. Die dabei aufgenommene Informa- 
tion wind helligkeitskodiert in einem Bild oder - in drei Di- 
mensionen — als Volumen dargestellt. Als bildgebende Sy- 65 
steme kdnnen unterschiedliche Vorrichtungen verwendet 
werden. Nebcn Vorrichtungen zur Erfassung histologischcr 
Bildsequenzen mit lichtmikroskopischer Information ist es 



auch moglich, makroskopische Information aus Rontgen- 
Computertomographie-(CT), Kemspin tomograph ie-(MRT), 
funktioncllc MRT-(fMRT), Positron-Emissionstomogra- 
phie-(PET) oder Single-Photon-Emission-Computer-Tomo- 
graphie-(SPECT) Bildsequenzen zu verwenden. 

Unter Verwendung dieser bildgebenden Verfahren bctrifft 
die Erfindung insbesondere die vollautomatische Berech- 
nung von Raumdifferenzen biologischer Objekte zur raum- 
lichen Normierung. Die Berechnung der Raumdifferenzen 
erfolgt hierbei flir jedes Volumenelement, mit einer Kombi- 
nation aus einem Bewegungsmodell, Mehrgitter- Verfahren, 
externer Krafte und einem Mechanismus zur Vergroberung 
der Volurnina. Die Gesamtheit der Raumdifferenzen (Defor- 
mationsfelder) wird von einem Speicher erfaBt, der mit einer 
Visualisierungseinheit verbunden ist. Die Deformationsfel- 
der konnen auch ftir die Integration der Bildinformation ver- 
schiedener, bildgebender Verfahren verwendet werden, d. h. 
fur die tJberlagcrung der mit vcrschiedenen Verfahren gc- 
wonnenen Bildinformationen in ein gemeinsames Koordi- 
natensystem. Durch das hier vorgestellte Verfahren wird es 
moglich, die Objekte mit hochstmoglicher geometrischer 
Genauigkeit miteinander zur Deckung zu bringen, ohne daB 
eine Zuordnung von entsprechenden Volumenelementen mit 
einer meBtechnischen Vorrichtung erforderlich wird. Dieses 
Verfahren ermittelt automatisch die Distanz fur jedes Volu- 
menelement des Ausgangsobjekts zum Referenzobjekt, was 
mit bestehenden MeB- oder Zuordnungsmethoden aufgrund 
der groBen Anzahl der Volumenelemente (mindestens ca. 8 
Millionen in der Medizintechnik) nicht moglich ist. Insbe- 
sondere ist es mit klassischen Verfahren zum Ausmessen 
dreidimensionaler Objekte nicht moglich, interne Volumen- 
elemente des Ausgangsobjekts in dem Referenzobjekt wie- 
derzuerkennen, weil sie i. A. entweder nicht zuganglich 
oder auch nicht besonders ausgebildet sind, so daB eine Zu- 
ordnung nicht moglich ist. Es ist bereits bekannt, Anpas- 
sungsmethoden bereitzustellen, die interaktiv mit wenigen 
Stutzstcllcn und cingeschrankter Zuordnungsmoglichkeit 
(Roland et al., 1994; Schiers et al., 1986; Bookstein, 1989; 
Rohr et al., 1996, vgl. auch Literaturverzeichnis am Ende 
der Schrift) oder auch automatisch (Friston et al., 1995; Lin- 
deberg, 1995; Collins et al., 1995) aber nur Referenzpunkt-, 
Oberflachen- oder Kontur-basierend (Davatzikos, 1996; 
Thompson et al. 1996) oder global (Bajcsy et al., 1989; Mil- 
ler et al., 1993; Bro-Nielsen et ab, 1996; Thirion, 1996) ar- 
beiten. Diese bekannten Verfahren haben die folgenden un- 
vermeidbaren Einschrankungen, die eine fur die Medizin- 
technik relevante vollautomatische und volumenclementgc- 
naue Anpassung unmoglich machen: 

1 . Manuell unterstiitzte Verfahren sind nicht reprodu- 
zierbar, so daB eine Wiederholung der Messung zu un- 
tcrschiedlichen Ergebnissen fuhren wiirdc. Diese Ver- 
fahren lassen sich sinnvoll nur bei eingeschrankter Zu- 
ordnung (z. B. globale Drehung) verwenden, weil an- 
dernfalls die Handhabung einer groBen Anzahl von Re- 
ferenzpunkten fur eine genauere Zuordnung praktisch 
nicht durchfuhrbar ist. 

2. Die Verwendung von nur Kontur- oder Oberflachen 
fuhrt - neben der ebenfalls notigen manuellen Unter- 
stiitzung - zu Interpolationsproblemen entsprechender 
Oberflachen oder Konturen zwischen Ausgangs- und 
Referenzobjekt. Zudem ist die Zuordnung i. A. nicht 
reproduzierbar. 

3. Es wird die Voraussetzung von nur kleinen Defor- 
mationen gemacht, wornit eine vollstandige Zuordnung 
nicht moglich wird. 

4. Unter der Voraussetzung von groBen und komplc- 
xen Raumunterschieden ist der numerische Rechenauf- 
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wand so gro8, daB die bisher bekannten Verfahren ent- 
weder auf eine vollstandige Zuordnung verzichten 
miissen oderes werden zur Minimierunng des Rechen- 
aufwandes Filter verwendet. Allerdings ist es mil Hilfe 
dieser Filter nicht in ausreichender Weise moglich, eine 5 
volumenelementgenaue Anpassung zu erreichen, weil 
diese filterbasierenden Verfahren entweder heuristisch 
arbeiten, womit die Bestimmung exakter Dcformati- 
onsfelder nicht mehr moglich wird oder die Umsetzung 
ist fur eine i. A. dreidimensionale Anwendung in der io 
Praxis wegen des groBen Rechenaufwandes immer 
noch nicht durchfiihrbar. Hierfur ist entscheidend, daB 
die Anzahl der aufzuwcndendcn Arbcitsschritte O die- 
ser filterbasierenden Verfahren im besten Fall mit 0(N 
logN) anwachst (N ist die Gesamtanzahl der Volumen- 15 
elemente). 

Die vorliegende Erfindung vermeidet die o. a. Nachteile 
(1-4) vollstandig. Bei dem in Schormann et al. (1996) be- 
schriebenen Verfahren handelt es sich um einen weiteren 20 
Vorschlag zur Bestimmung von relativ zu der vorliegenden 
Erfindung mit einem Bewegungsmodell - nur kleinen geo- 
metrischen Raumunterschieden. Zudem beriicksichtigt die- 
ser Vorschlag keinen Mechanismus zur Vergroberung der 
Objekte zur weiteren Verbesserung der Verarbcitungsge- 25 
schwindigkeit, keine modifizierten, extemen Krafte zur 
Steigerung der Anpassungsqualitat. Zusatzlich wird in die- 
ser Arbeit und in Henn et al. (1997) ein nicht-lineares Mo- 
dell beschrieben, das gegenuber dem linearen Modell Nach- 
teile hinsichtlich der numerischen Konvergenz, numerischer 30 
Instability und Konsistenz aufweist (s. u.). In diesem Zu- 
sammenhang ist zu erwahnen, daB es keine "Standard- 
Mehrgitter- Verfahren" fur samtliche, in diesem Sachzusam- 
menhang auftretenden Differentialgleichungen existieren. 
Sondern die Mehrgitter- Verfahren bcstehen aus vielen ein- 35 
zelnen Komponenten, die auf die jeweilige Problemstellung 
zugeschnittcn werden miissen. Der Unterschicd fur die L6- 
sung eines linearen bzw. nicht- linearen Modells ist daher 
durch die unterschiedlichen Anforderungen der Modelle 
selbsl vorgcgcbcn und beinhaltet daher eine eigene, an das 40 
Modell angepaBte Struktur (Art und Anzahl der Post- und 
Prarelaxation, Restriktion, Interpolation, V-, W-, F-Cycle, 
Gittertiefe, vgl. Press WH, Flannery BP, Teukolsky SA, Vet- 
terling WT: Numerical recipes in C. Cambridge: Cambridge 
U. P. (1988), S. 871 ff), womit es erforderlich war, fur das li- 45 
neare System ein neues Mehrgitterverfahren zu entwickeln. 

Daruber hinaus unterscheidet sich das Bewegungsmodell 
von dem in Christensen et al. ((1997), S. 868, Abschnitt 3) 
vorgeschlagenen Fliissigkeitsmodell in den folgenden Punk- 
ten: 50 

a) In der Arbeit von Christensen et al. ((1997), (1996)) 
wird die Umsetzung zu einem Fliissigkeitsmodell 
durch zusatzliche Losung einer gewohnlichen Diffe- 
rentialgleichung (GL 10 in Christensen et al. ((1997)) 55 
erreicht. Wie aus den Gleichungen 10 und 12 dieser Ar- 
beit crsichtlich ist, muB fur jeden Zcitschritt die Ablei- 
tung fur jede Raumrichtung (der Gradient des Defor- 
mationsfeldes Vu Vu) berechnet werden. Diese Be- 
rechnung wird durch das von uns vorgestellte Bewe- 60 
gungs modell vollstandig vermieden, was unter der Be- 
riicksichtigung der groBen Datenmengen und insbeson- 
dere der Zielsetzung einer schnellen Berechnung der 
Raumunterschiede einen erheblichen Rechenzeitvor- 
teil bietet. 65 

b) Zudem unterscheidet sich das Bewegungsmodell 
grundsatzlich von den in Christensen et al. (1996), 
(1997); Schormann et al. (1996); Henn et al. (1997) 
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und Haller et al. (1996) beschriebenen Ausfuhrungen: 
(i) in der linearisierten Ausfiihrung des Modells: 

U8u (n >) = T(U gc / D >) 

wobci L der Differentialoperator L = uA + (u + 
X) V V • , u, X die Lam^'schen Konstanten,T die grau- 
wertbasierende, externe Kraft, 8u (n) die n-te lokale Ver- 
schiebung im Modell und u gcs (n) die bereits berechnete, 
gesamte Bewegung eines beliebig herausgegriffenen 
Volumenelcmentes sind. Die oben beschricbene Linea- 
risierung des Modells durch Einfuhrung zeitunabhan- 
giger, infinitesimaler Wegstrecken 5u ( ^, die in ihrer 
Gesamtheit die Bewegung der Volumenelemente be- 
schreiben, stellt einen gravierenden Unterschied zu den 
in den o. a. Arbeiten genannten Verfahren da (vgl. auch 
Fig. 1). Bei dem in Christensen et al. bevorzugten zeit- 
abhangigen Modell muB, im Vergleich hierzu, eine zu- 
satzliche gewohnliche Differentialgleichung gelost 
werden. 

Das in Schormann et al. (1996) und Herrn et al. (1997) be- 
schricbene Verfahren beinhaltet ein nicht-lineares Modell, 
das sowohl im Differentialoperator als auch bei der externen 
Kraft T = T(u) Abhangigkeiten in den Raumunterschieden u 
aufweist: 

L0D=T(IT) 

wobei eine Taylorreihen-Entwicklung von T(u) auf die 
Nichtlinearitat bzgl. u fuhrt. Diese nicht-linearen Modelle 
weisen folgende Nachteile auf: Zum einen sind sie sehr an- 
fallig gcgeniibcr numerischen Instabilitaten, die praktisch 
unvermeidbar durch Rauschen der bildgebenden Systeme 
verursacht werden. Zum anderen stellt die sehr geringe nu- 
merische Konvergenz bei diesen nicht-linearen Systemen 
ein Problem dar und zudem ist die Konsistenz, d. h., die 
Korrelation homologer Strukturen im Ausgangs- und Ziel- 
volumen problematisch. 

Das Bewegungsmodell unterscheidet sich daher grund- 
satzlich von dem Stand der Technik hinsichtlich der Zeitun- 
abhangigkeit und der Linearisierung. 

Zudem ist es mit Hilfe der bekannten Verfahren nicht in 
ausreichender Genauigkeit moglich, bei beliebigen Objek- 
ten automatisch, d. h. ohnc manuelle Unterstiitzung, Anpas- 
sungsergebnisse zu erzeugen, bei denen jedes Volumenele- 
ment des Ausgangsobjekts einem entsprechenden Volumen- 
element des Referenzobjekts zugewiesen wird. 

Es ist daher Aufgabe der vorliegenden Erfindung, ein Ver- 
fahren der einleitend genannten Art zu verbesscrn, so daB 
der extrem hohe numerische Rechenaufwand bei gleichblei- 
bend hoher Anpassungsgiite durch ein linearisiertes Bewe- 
gungsmodell und durch ein auf dieses Modell angepaBtes 
Mehrgitter- Verfahren erheblich reduziert wird, womit die 
Berechnung von Deformationsfeldcrn ohnc intcraktive 
Hilfe ermoglicht wird. Die Berechnungszeiten in drei Di- 
mensionen fiir das Fliissigkeitsmodell werden in der Litera- 
tur (vgl. Christensen et al., IEEE Computer, 29(1): 32-38, 
1996) mit ca. 7 Tagen angegeben, wohingegen das hier vor- 
gestellte Verfahren - in einer parallelisierten Ausfiihrung - 
ca. eine halbe Stunde benotigt. 

Diese Aufgabe wird erfindungsgemaB durch ein Verfah- 
ren mit den Merkmalen des Anspruchs 1.) dadurch gelost, 
daB zur Berechnung der Deformationsfelder die Objekte zu- 
nachst von bildgebenden Verfahren in einem computcr-ge- 
stutzten System erfafit werden. Nach der Erfassung wird fur 
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jedes Volurnenelement des Ausgangsvolumens die Bewe- 
gung auf ein entsprechendes Volurnenelement des Referenz- 
volumcns bestimmt. Hierfur werden die Objekte zunachst 
zum elastischen Stoff gewandelt, indem jedes Volurnenele- 
ment uber Fedem mit seinen nachsten Nachbarn verbunden 5 
wird. Dicser Vorgang enlsprichl einer Diskrctisierung des 
ursprunglich kontinuierlichen Problems. Durch Anbringen 
einer auBcren Kraft, die gerade erforderlich ist, urn die 
Raumdifferenz zu einem Referenzvolumenelement zu iiber- 
winden, deformiert sich dieses aus Volumenelementen be- 10 
stehende Gitter, wobei allerdings die enlstehende Deforma- 
tion nur ein KompromiB zwischen der exteraen Kraft und 
den aus den Fedem resultiercnden, internen Kraften ist (Ela- 
stizitatsmodell). Zur Berechnung der vollstandigen Raum- 
differenz wird daher an der neuen Zwischenposition eine 15 
neu berechnete, exteme Kraft aufgewendet bis schlieBlich 
der Endzustand, eine vollstandige Deckung mit dem ent- 
sprechenden Referenzvolumenelement, erreicht ist. 

Durch diese sukzessive, d. h. linearisierte Anwendung der 
Krafte, kann die gesamte Bewegung eines beliebigen Volu- 20 
menelementes bestimmt werden, womit groBe Deformatio- 
nen berechnet bzw. korrigiert werden konnen. Die Gesamt- 
bewegung der Volumenelemente wird einer bevorzugten 
Ausfuhrungsform erfindungsgemaB durch sukzessive An- 
wendung neuer, externer Krafte erreicht, die sich aus den ge- 25 
anderten Raumdifferenzen des vorhergehenden Ortsinter- 
valls ergeben. Sie kann aber auch durch Losung einer zeitab- 
hangigen Bewegungsgleichung fur jedes Volurnenelement 
erreicht werden, was jedoch numerisch aufwendiger ist 
(Christenscn ct ah, 1994). Eine weitere Aufgabe der vorlie- 30 
genden Erfindung ist es, geeignete Krafte zwischen der 
Form des Ausgangs- und Referenzvolumens zu bestimmen, 
mit denen der zuruckgelegte Weg jedes Volumenelementes 
des Ausgangsvolumens auf das Referenzvolumen gesteuert 
werden kann. Diese Aufgabe wird erfindungsgemaB da- 35 
durch gelost, daB die quadrierte Grauwertdifferenz zwi- 
schen dem Referenzvolumenelement und dem um eine 
kleine Raumdifferenz verschobenen Ausgangsvolumenele- 
ment als auch die quadrierte Grauwertdifferenz zwischen 
Ausgangsvolumenelement und dem um eine kleine Raum- 40 
differenz verschobenen Referenzvolumenelement minimiert 
wird. Hiermit wird aus der Grauwertdifferenz gerade die 
Starke, der Ort und die Richtung der Kraft bestimmt, die fur 
jedes Volurnenelement die sich einstellende Bewegungs- 
bahn eindeutig beschreibt. Durch die Symmetrisierung zwi- 45 
schen Ausgangs- und Referenzvolumen wird die Kraft in ei- 
ncn Antcil aufgespalten, der das Ausgangsvolumen in das 
Referenzvolumen driickt und einen Anteil der das Aus- 
gangsvolumen in das Referenzvolumen zieht. Eine andere 
denkbare Alternative ist eine Kraft, die nur aus dem 1. An- 50 
teil der symmetrischen Kraft besteht, wobei aber eine uner- 
wiinschtc Vcrschlechterung des numerischen Aufwandes 
und der Anpassungsgenauigkeit erreicht wird. Altemativ 
zur symmetrisierten Kraft werden die Krafte erfindungsge- 
maB unempfindlich gegeniiber Rauschen bildgebender Sy- 55 
steme gemacht, indem die Starke der Krafte zunachst qua- 
dratisch mit den Grauwertdiffcrerizen anwachsen und fur 
groBere Differenzen linear zunehmen. 

Eine weitere Aufgabe der vorliegenden Erfindung ist es, 
die Deformationsfelder mit moglichst geringem numeri- 60 
schen Aufwand ohne Verlust an Genauigkeit zu berechnen. 
Diese Aufgabe wird erfindungsgemaB mit. Hilfe von finiten 
Differenzen oder finiten Elementen als Diskretisierer und 
auf das Modell angepaBten numerischen Mehrgitter- Verfah- 
ren, dem CS (Correction Scheme) oder dem FAS (Full Ap- 65 
proximation Scheme) gelost. Hiermit wird es ermoglicht, 
daB der numcrischc Aufwand nur linear mit der Anzahl der 
Volumenelemente anwachst (0(N)). Die bisher besten be- 



kannten globalen Filter- Verfahren wachsen mit 0(N logN) 
(Fliissigkeitsmodell Christensen et al., 1997: 0(N 2 )) an, wo- 
mit bei der groBcn Anzahl der Volumenelemente (klinische 
Anwendungen: ca. 8 Millionen) ein erheblicher Zeitnachteil 
erwachst, der diese Verfahren fUr die praktische Anwendung 
mil bisher bekannten Rcchncrleistungcn unbrauchbar 
macht. Bei der vorliegenden Erfindung wird das Ergebnis 
nicht approximicrt, sondern kann mit einer gewahlten nu- 
merischen Genauigkeit bestimmt werden. Als Ergebnis wird 
der von jedem Volurnenelement zuruckgelegte Weg in ei- 
nem MeBwertspeicher erfaBt. Daran anschlieBend wird die 
punktweise Verschiebung von dem Ausgangsvolumen mit 
Hilfe des Deformationsfcldes durchgefuhrt und das Ergeb- 
nis von einem Speicher erfaBt oder einer Visualisierungsein- 
heit zugefuhrt. 

Eine weitere Aufgabe der Erfindung ist es, die automati- 
sche Zuordnung globaler, geometrischer Strukturen inner- 
halb der Objekte zu verbessem und cine weitere Steigerung 
der Verarbeitungsgeschwindigkeit zu erreichen. Diese Auf- 
gabe wird erfindungsgemaB mit Hilfe von einem Mechanis- 
mus zur Vergroberung der Volumina erreicht. Bei diesem 
Verfahren werden mehrere benachbarte Volumenelemente 
zu einem neuen Volurnenelement zusammengefaBt, z. B. 
starr verbunden. Dieses Verfahren wird mehrfach wieder- 
holt, so daB cine Hierachie von grober werdenden Volumina 
entsteht. Das Deformations feld wird zunachst fur das grob- 
ste Volumen bestimmt. Die Deformationsfelder werden auf 
das nachst feinere Volumen interpoliert bis schlieBlich die 
originale Auflosung erreicht ist. Durch die Anwendung auf 
zunachst grobe Volumina gelangt die lokale Zuordnung von 
Volumenelementen in den Bereich globaler Raumdifferen- 
zen, so daB hiermit die Anpassung grober Strukturen und 
damit uber groBere Raumdistanzen erfolgen kann. Durch 
diese Reduktion der Gesamtanzahl der Volumenelemente 
wird eine Steigerung der Verarbeitungsgeschwindigkeit er- 
reicht. Gleichzeitig wird die automatische Zuordnung geo- 
metrischer Strukturen insbesondcre bei raumlich komplcxen 
Strukturen der Objekte verbessert, da die Zuordnung von 
groben nach feinen Strukturen systematisch aufgebaut wird. 
Mit der vollautomatischen Berechnung von Deformations- 
feldem steht eine Erfindung zur Verfugung, die es im Be- 
reich der Medizintechnik moglich macht, eine prazise Loka- 
lisation der Daten verschiedener bildgebender Verfahren 
vorzunehmen. Es stellt ein neues, bisher nicht verfugbares 
Werkzeug insbesondere zur Analyse der Cortexorganisation 
beim Menschen dar. Eine solche Analyse ist nicht moglich, 
wenn nicht vorher das AusmaB der durch normale interindi- 
viduelle Variabilitat bedingten Strukturvariation bekannt ist. 
Diese geometrische Normalitat von Strukturen innerhalb der 
Objekte kann mit Deformationsfeldern berechnet werden. 
Hierfur werden entweder die Strukturen einer Population 
von Objektcn mit Hilfe der Deformationsfelder in ein ge- 
meinsames Koordinatensystem transformiert, wobei dieses 
Verfahren eine Analyse der Uberlagerungswahrscheinlich- 
keiten dieser Strukturen ermoglicht. Oder aus den Deforma- 
tionsfelder einer Population von Objekten lassen sich GauB- 
Felder berechnen. Hiermit ist es uber einen Chi-Quadrat- 
Test bezuglich des Mahanalobis-Abstandes, der sich aus El- 
lipsoiden gleicher Wahrscheinlichkeitsdichten fiir jedes Vo- 
lurnenelement berechnet, moglich, Signifikanzniveaus zur 
Normalitat von Objekten anzugeben. Die Durchfuhrung die- 
ser Analyse soil in den Aufbau einer speziellcn Datcnbank 
miinden, in der alle Daten in einem sogenannten "menschli- 
chen Referenzgehirn" zusammengefaBt sind. 

Unter Beriicksichtigung individueller anatomischer Un- 
terschiede lassen sich unter Zuhilfenahme dieser Datenbank 
pathologische Verandcrungcn quantitativ und qualitativ cr- 
fassen und adaquat behandeln. In der klinischen Diagnostik 
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konnen damit die Grenzen der Normalitat aktueller Patien- 
tendaten bzgl. eines Referenzobjektes validiert werden, was 
z. B. von groBcr Bcdcutung fur stcreotaktischc Opcrationen 
ist. 

5 

Bcschreibung der Zeichnungen 

Die folgenden Zeichnungen dienen dem bcsseren Ver- 
standnis der Erfindung, wobei weitere Einzelheiten den 
nachfolgenden, ausflihrlichen Beschreibungen zu entneh- 10 
men sind. In den Zeichnungen zeigt Fig. 1 die Bewegung 
der "wahren" Bahn eines Volumenelementes, der Approxi- 
mationcn durch das lineare Bewegungsmodell 5u*i und ci- 
nem nicht-linearen Modell (geschwungene Kurve). Fig. 2 
ist eine schematische Zeichnung einer Ausfuhrungsform der 15 
vorliegenden Erfindung. Fig. 3 zeigt eine zweidirnensionale 
Anwendung des Deformationsfeldes zur Verkniipfung von 
MR- (a) und cincm entsprechenden histologischen Scrien- 
schnitt (b) zur Integration multimodaler Bildinformation. 
Erkennbar ist, daB trotz unterschiedlicher Schnittrichtungen 20 
bei der Datenaufnahme und unvermeidbarer Herstellungsar- 
tefakte bei der histologischen Aufbereitung eine volumen- 
elementgenaue Zuordnung erreicht ist. Fig. 4 zeigt eine 
zweidirnensionale Anwendung des Deformationsfeldes zur 
Verkniipfung von zwei verschiedenen Gehirnen (a, b), die 25 
mit Hilfe der MR-Tomographie aufgenommen wurden. Um 
die raumlichen Unterschiede in beiden Gehirnen vor (a) und 
nach (c) der Anpassung sichtbar zu machen, wurde die au- 
Bere Kontur des Referenzschnittes (b) in (a) und (c) iiberla- 
gert. Fig. 5 zeigt das Deformationsfeld (a), das gcrade die in 30 
Fig. 4 gezeigte Ausgangsstruktur auf die Referenzstruktur 
deformiert. Um geometrische Veranderungen auch innerer 
Strukturen sichtbar machen zu konnen, wird das Deformati- 
onsfeld auf ein regelmaBiges Gitter angewendet, wobei Be- 
reiche des Gehirns hervorlreten, in denen es expandiert bzw. 35 
komprimiert wird. Darunter sind die Krafte fiir die x- (b) 
und y-Richtung (c) gczeigt, wobei die Starke und Richtung 
in Grauwerte so kodiert ist, daB in (b) die Kraft, die nach 
links zeigt in weiB (schwarz nach rechts) und in (c) die 
Kraft, die nach oben zeigt in weiB (schwarz nach unten) dar- 40 
gestellt ist. Der konstante Grauwert am Rand des Bildes 
markiert eine Kraft der Starke 0. Die resultierende Kraft in 
jedem Bildpunkt setzt sich aus der vektoriellen Summe die- 
ser Krafte zusammen. 

Fig. 6 zeigt ein weiteres Anwendungsbeispiel (Ausgangs- 45 
objekt vor (a) und nach (c) der Anpassung, (b) Referenzob- 
jekt), bei dem besonders starke Deformationen im rechten 
Bereich des Objekts auftreten (vgl. Deformationsfeld (d)). 
Fig. 7 zeigt eine dreidimensionale Anwendung des Verfah- 
rens an zwei individuellen Gehirnen. In der Ausgangssitua- 50 
tion ist die sehr unterschiedliche Form der Oberflache (Sul- 
cusmuster) erkennbar. Nach der Deformation sind die Volu- 
mina beider Gehirne, insbesondere das Sulcusmuster iden- 
tisch. Die Erfassung der geometrischen Ausdehnung eines 
dreidimensionalen Objekts kann nach Fig. 2 mit Hilfe von 55 
bildgebenden Verfahren der Medizin (2) erfolgen, die mit 
mindestens einem Dctcktor ((1), z. B. Spule, CCD) ausge- 
stattet sind. Die dabei erfaBten Signale werden einem A/D 
Wandler (3) zugefuhrt und die digital gewandelten Signale 
des Volumens 1 und des Volumens 2 in einem Speicher (4a), 60 
(4b) gehalten. Bei der Erfassung der Objekte mit bildgeben- 
den Verfahren sind die Volumina mit Hilfe von stercotakti- 
schen Methoden so zu positionieren, daB eine globale ent- 
sprechende Ausrichtung der Volumina gewahrleistet ist. Bei 
fehlenden stereotaktischen Moglichkeiten wird die Position 65 
und Ausrichtung global durch Ausfuhrung einer affinen Be- 
wegung des Ausgangsobjckts korrigiert (5). Hierbei ge- 
wahrleistet eine affine Bewegung i. A. eine hohere Anpas- 



sungsgiite als die Moglichkeit der stereotaktischen Fixie- 
rung, weil neben der Position auch noch zusatzlich die 
GroBc korrigiert wird. Als nachstcr Schritt erfolgt cine 
Wandlung der Objekte als elastisches Material (6), indem je- 
des Volumenelement uber Federn mit seinen nachsten Nach- 
barn verbunden wird. Hierdurch wird das Volumen in ein 
aus Volumenelementen bestehendes Gitter gewandelt. Jedes 
Volumenelement bcsteht aus einem bestimmten Grauwert. 
Die Krafte, die zur Uberwindung der RaumdifTerenzen er- 
forderlich sind, ergeben sich - nach Anwendung von Extre- 
malprinzipien der theoretischen Physik (Prinzip von Hamil- 
ton) - aus den Grauwertdifferenzen beider Volumina fur je- 
des Volumenelement, was den Vorteil hat, daB keinc ent- 
sprechenden PaBpunkte in beiden Volumina bestimmt wer- 
den m iis sen. Sondern vielmehr werden Volumenelemente 
des Ausgangsvolumens so verschoben, daB sich die Grau- 
wertdifferenz zwischen beiden Volumina vermindert. Bei 
der elastischcn Wandlung (6) erfolgt diescr Vorgang cinmal, 
wobei die resultierende Verschiebung der einzelnen Volu- 
menelemente ein KompromiB ist zwischen der internen 
Kraft, die sich aus dem Federsystem berechnet und der ex- 
ternen Kraft, die sich aus der Grauwertdifferenz ergibt. In ei- 
nem nachsten Schritt werden mehrere Volumenlemente zu 
einem neuen Volumenelement zusammengefafit (7), indem 
die Fcderverbindung zwischen ihnen durch eine starre Vcr- 
bindung ersetzt wird und der neue Grauwert sich z. B. aus 
dem Mittelwert der zusammenzufassenden Volumenele- 
mente ergibt. Dieses Verfahren wird mehrfach wiederholt, 
so daB eine Hierarchie aus grober werdenden Volumina ent- 
stcht. Das Deformationsfeld (9) wird zunachst auf dem 
grobsten Volumen mit Hilfe von Mehrgitter- Verfahren (8) 
bestimmt. 

Eine weitere Verbesserung der Anpassungsgute erfolgt 
durch die Berechnung der Gesamtbewegung (10) durch suk- 
zessive Anwendung des Deformationsfeldes (9) undLosung 
mit dem Mehrgitter- Verfahren (8). Hierdurch en ts tent eine 
Sequenz von Volumina, die dem Referenzvolumen iramcr 
mehr ahnelt, bis sie in dieser Auflosungsstufe identisch sind. 
Das resultierende Deformationsfeld wird auf das nachst fei- 
nere Volumen inlerpoliert, der ProzcB (8)-(9)-(10) auf jeder 
Auflosungsstufe wiederholt, bis schlieBlich die originale 
Auflosung erreicht ist. Durch die Anwendung auf zunachst 
grobe Volumina gelangt die lokale Zuordnung von Volu- 
menelementen in den Bereich globaler RaumdifTerenzen, so 
daB hiermit die Anpassung grober Strukturen und damit 
uber groBere Raumdistanzen erfolgen kann, womit auch 
gleichzeitig die automatische Zuordnung raumlich komplc- 
xer Strukturen verbessert wird, da die Zuordnung von gro- 
ben nach feinen Strukturen systematise h aufgebaut wird, bis 
schlieBlich der Endzu stand, eine vollstandige Obereinstim- 
mung mit dem entsprechenden Referenzvolumen erreicht 
ist. Als bildgebende Systeme (2) konnen unterschiedliche 
Vorrichtungen verwendet werden. Neben Vorrichtungen zur 
Erfassung histologischer Bildsequenzen mit lichtmikrosko- 
pischer Informadon ist es auch moglich, makroskopische 
Information aus Rontgen-Computertomographie-(CT), 
Kcrnspintomographic-(MRT), funktionelle MRT-(tMRT), 
Positron-Emissions-Tomographie-(PET) oder Single-Pho- 
ton-Emission-Computer-Tomographie- (SPECIE B i Idse- 
quenzen zu verwenden. Hierbei kommt der Integration mul- 
timodaler Bildinformation aus histologischen und MR-Bild- 
sequenzen des idendschen Objekts eine besondcrc Bedcu- 
tung zu, weil hiermit die nur makroskopische Information 
der MR-Tomographie durch die lichtmikroskopische Infor- 
mation erganzt werden kann, sofern das histologische Volu- 
men mit dem entsprechenden MR- Volumen durch die vor- 
liegende Erfindung vollstandig zur Dcckung gcbracht ist 
(Fig. 3). Hiermit wird die makroskopische Auflosung der 
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MR-Tbmographie um einen Faktor 1000 verbesserl, womit 
die Abgrenzung mikroskopisch definierter, anatomischer 
Areale ermoglicht wird. Die Erfindung ist dahingehend kon- 
zipiert (Fig. 2), daB es nicht erforderlich wird, fur jede 
Struktur ein Deformationsfeld neu zu berechnen, sondern 5 
jede separate Struktur 1 . . . N des Ausgangsobjcktes kann 
mit Hilfe des einmal berechneten Deformationsfeldes auf 
die Refercnz deformicrt wcrden. 

Aus der Darstellung Fig. 4 ist ersichtlich, daB die sehr un- 
terschiedliche Morphologie verschiedener Gehirne durch 10 
die Erfindung vollautomatisch korrigiert werden kann. Die 
dazu erforderlichen Krafte sind in Fig. 5 dargestellt, womit 
vcrdeutlicht wcrden soli, daB fur jedes Volumcnelement eine 
Kraft einer bestimmten Starke und Richtung berechnet wird. 
AuBerdem ist auch ersichtlich, daB die Kraft erfindungsge- 15 
maB aus zwei Anteilen besteht: einem Anteil, der das Aus- 
gangsvolumen in das Referenzvolumen driickt und einem 
Anteil der das Ausgangsvolumen in das Referenzvolumen 
zieht. Dies ist erkennbar an der ausgebildeten Uberlagerung 
der beiden Objekte fur die x- und y-Richtung. Fig. 5 zeigt 20 
das visualisierte Deformationsfeld, das auf ein regelmaBiges 
Gitter angewendet wurde, wobei veranschaulicht wird, daB 
groBe und komplexe Deformationen auch irn Innern des Ob- 
jekts vorhanden sind (Fig. 6), was an den teilweise stark ex- 
pandierten odcr komprimicrten urspriinglich quadratischen 25 
Kastchen deutlich wird. Fig. 7 zeigt eine dreidimensionale 
Realisierung der Erfindung. 
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Patentanspriiche 

1. Verfahren zur Quantifizierung und Darstellung von 
2- oder 3-dimensionalen Raumunterschieden fur die 
medizinische Diagnostik, Forschung oder Medizin- 
technik, mit dem die relativen Positionen der aus der 
Digitalisierung bildgebender Verfahren erfaBten Volu- 
menelemente mindestens eines Objekts zu mindestens 
einer Referenzstruktur und i. A. zu samllichen Volu- 
menelementen dieser Referenz berechnet, gespeichert 
und dargestellt werden, dadurch gekennzeichnet, daB 
die aus der Digitalisierung bildgebender Verfahren re- 
sultierenden Volumenelementc mindestens zweier Ob- 
jekte von einem computer- gestiitzten System erfaBt 
werden, wobei die zuriickgelegten Wege der Volumen- 
elementc des Ausgangsvolumens auf die entsprechen- 
den Referenz- Volumenelemente durch ein linearisier- 
tes Bewegungsmodell, ein auf dieses Bewegungs mo- 
dell angepaBtes numerisches Mehrgitter- Verfahren, 
grauwertbasierenden Kraften und/oder mit mehreren 
Auflosungsstufen der dig i tali sierten Volumina berech- 
net werden. 

2. Verfahren nach Anspruch 1, dadurch gekennzeich- 
net, daB der von jedem Volumenelement zuriickgelegte 
Weg (Deformationsfeld) in einem MeBwertspeicher er- 
faBt wird. 

3. Verfahren nach den Anspriichen 1 und 2, dadurch 
gekennzeichnet, daB die Verschiebung der Volumen- 
elemente von dem gesamten oder Teilen des Ausgangs- 
volumens mit Hilfe des Deformationsfeldes durchge- 
fuhrt wird. 

4. Verfahren nach einem der Anspriiche 1 bis 2 da- 
durch gekennzeichnet, daB cine Visualisierung des De- 
formationsfeldes mit Hilfe von einem regelmaBigen 
Gitter einer Visualisierungseinheit (Druckermedium, 
Bildschirm) zugefuhrt wird. 

5. Verfahren nach einem der Anspriiche 1 bis 3, da- 
durch gekennzeichnet, daB die Verschiebung der Volu- 
menelemente von dem gesamten oder Teilen des Aus- 
gangsvolumens einer Visualisierungseinheit (Drucker- 
medium, Bildschirm) zugefuhrt wird. 

6. Verfahren nach einem der Anspriiche 1 bis 5, da- 
durch gekennzeichnet, daB die Wahrschcinlichkcit des 
Auftretens von Strukturen einer Population von Objek- 
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ten durch Anwendung der Deformalionsfelder auf die 
Strukturen selbst oder die Wahrscheinlichkeiten durch 
Verwcndung der Deformation sfelder bcrechnet wer- 
den. 

7. Verfahren nach einem der Anspriiche 1 bis 6, da- 5 
durch gekennzeichnet, daB ein Elastizitatsmodell ver- 
wendet wird, bei dem das resultierende Deformations- 
fcld einmal und nicht sukzessiv berechnet wird. 

8. Verfahren nach einem der Anspriiche 1 bis 7, da- 
durch gekennzeichnet, daB die Krafte sich aus einer 10 
symmetrisierten Grauwertdifferenz zwischen Aus- 
gangs- und Referenzvolumen ergeben. 

9. Verfahren nach einem der Anspriiche 1 bis 7, da- 
durch gekennzeichnet, daB die Starke der Krafte zu- 
nachst quadratisch mit der Grauwertdifferenz anwachst 15 
und fur groBere Differenzen Linear zunimmt. 

10. Verfahren nach einem der Anspriiche 1 bis 9, da- 
durch gekennzeichnet, daB das Mehrgitter- Verfahren 
als Correction Scheme (CS) oder als Full Approxima- 
tion Scheme ausgelegt ist, um das linearisierte Bewe- 20 
gungsmodell zu losen. 

11. Verfahren nach einem der Anspriiche 1 bis 10, da- 
durch gekennzeichnet, daB die "Oberfuhrung (Diskreti- 
sierung) in ein Elastizitats- oder Bewegungsmodells 
mit Hilfe fin iter DifTcrcnzen oder finiter Elemente vor- 25 
genommen wird. 

12. Verfahren nach einem der Anspriiche 1 bis 11, da- 
durch gekennzeichnet, daB die Losung mit den Mehr- 
gitter- Verfahren und/oder das Verfahren zum Vergro- 
bern der Volumina auf Parallelrechnern bei gleichzeiti- 30 
ger Verwendung von mehreren Prozessoren durchge- 
fiihrt wird. 

13. Verfahren nach einem der Anspriiche 1 bis 12, da- 
durch gekennzeichnet, daB die Berechnung der Defor- 
mationsfelder fur verschiedene bildgebende Modalita- 35 
ten eines oder verschiedener Objekte erfolgt. 

14. Verfahren nach einem der Anspriiche 1 bis 13, da- 
durch gekennzeichnet, daB mindestens ein Objekt mit 
Hilfe der Computer-Tomographie und/oder der Ma- 
gnet- Resonanz- Tomographic (MRT) und/oder der 40 
funktionellen MRT und/oder von histologischen Seri- 
enschnitten und/oder der Positronen-Emissions- Tomo- 
graphic (PET) und/oder der Single-Photon-Emission- 
Computer- Tomographic (SPECT) und/oder der Ma- 
gnetoenzephalographie und/oder der konventionellen 45 
Rontgentechnik und/oder der Ultraschalltechnik meB- 
technisch erfaBt ist. 
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